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ABSTRACT 

We present a new method for probing the physical conditions and metal enrichment of the In- 
tergalactic Medium: the composite spectrum of Lya forest absorbers. We apply this technique to a 
sample of 9480 Lya absorbers with redshift 2 < z < 3.5 identified in the spectra of 13,279 high-redshift 
quasars from the Sloan Digital Sky Survey (SDSS) Fifth Data Release (DR5). Absorbers are selected 
as local minima in the spectra with 2.4 < TL yQ < 4.0; at SDSS resolution (s» 150 kms^ 1 FWHM), 
these absorbers are blends of systems that are individually weaker. In the stacked spectra we detect 
seven Lyman-series lines and metal lines of O VI, N V, C IV, C III, Si IV, C II, Al II, Si II, Fe II, Mg II, 
and O I. Many of these lines have peak optical depths of < 0.02, but they are nonetheless detected 
at high statistical significance. Modeling the Lyman-series measurements implies that our selected 
systems have total H I column densities JVhi ~ 10 154 cm~ 2 . Assuming typical physical conditions 
p/p = 10, T = 10 4 — 10 45 K, and [Fe/H]= —2 yields reasonable agreement with the line strengths of 
high-ionization species, but it underpredicts the low-ionization species by two orders of magnitude or 
more. This discrepancy suggests that the low ionization lines arise in dense, cool, metal-rich clumps, 
present in some absorption systems. 

Subject headings: galaxies: formation — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

The Lyman-a forest is a tracer of the diffuse mat- 
ter between galaxies known as the Intergalactic Medium 
(IGM), and at high redshift (2 < z < 4) this diffuse 
medium contains roughly 70% of the baryonic mass and 
occupies 90% or more of the volume. Galaxies form from 
the collapse of dark matter structures and corresponding 
accretion of gas from the IGM. But this is not a one- 
way process, and feedback is vital for a complete picture. 
Galactic outflows modify the kinetic and thermal energy 
of the medium and distribute the by-product of star for- 
mation: metals (e.g. Mac Low & Ferrara 1999; Madau 
et al. 2001; Oppenheimer & Dave 2006; Pieri et al. 2007; 
Pieri & Martel 2007). 

Metals have been observed in the high redshift Lya 
forest. In particular C IV (e.g. Meyer & York 1987 
Cowie et al. 1995; Songaila & Cowie 1996; Ellison et al 
2000; Schaye et al. 2003; Pieri et al. 2006) and O VI (e.g 
Schaye et al. 2000; Pieri & Haehnelt 2004; Simcoe et al 
2004; Bergeron & Herbert-Fort 2005; Aguirre et al. 2008; 
Frank et al. 2010; Pieri et al. 2010) provide prominent ab- 
sorbers and are commonly measured. Other species seen 
are Si IV, Si III, C III and N V (e.g. Schaye et al. 2003; 
Aguirre et al. 2004; Fechner & Richter 2009). These mea- 
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surements provide some indications of the regions of the 
Universe touched by mechanical feedback, but they can 
also provide a useful probe of the ionization properties of 
the medium (and so the physical conditions of the gas), 
and the abundance pattern. The main obstacle to sig- 
nificant progress on the latter two goals is a lack of large 
numbers of metal species to measure both from the same 
element (to measure ionization characteristics) and from 
different elements (to measure abundance patterns). 

In this paper, we introduce a new technique for the 
detection of weak metal lines — the composite spectrum 
of Lya forest absorbers. We apply this approach to the 
largest Lya forest dataset available: the Sloan Digital 
Sky Survey (SDSS) sample of QSO spectra. Pieri et al. 
(2010) demonstrated that the SDSS sample can provide 
precision measurements of weak lines by searching for 
O VI absorption associated with the bulk of the Lya for- 
est. Here we produce a blind search for absorption cor- 
related with the strong Lya forest absorbers; we identify 
and fit 19 metal lines. 

2. PRODUCTION OF COMPOSITE SPECTRA 

2.1. The Sample and Lyman-a Forest Line Selection 

We use the Sloan Digital Sky Survey Data Release 
5 (SDSS DR5) (Adelman-McCarthy et al. 2007), which 
provides 13,279 QSOs with useful forest coverage in the 
redshift range 2 < z < 3.5. The spectra are taken from 
the QSO Absorption Line Sample (QSOALS; York et al. 
2006) in the same manner described in Pieri et al. (2010) 
employed for DR3 spectra. The spectral resolution is 
wavelength dependent and varies from R — 1800 — 2200. 
In the following work we assume that R — 2000 at all 
wavelengths, and the error introduced by this approxi- 
mation is smaller than the quoted errors. 

We define "absorbers" in a simple way, selecting pixels 
that are lower in flux than their two neighboring pix- 
els and have optical depth 2.4 < tl vq < 4. We require 
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that our Lya absorbers are redward of the Ly/3 forest 
and blueward of 5000 kms -1 from the quasar Lya emis- 
sion redshift. In the stacked spectrum discussed in this 
paper we use only Lya absorbers in the redshift range 
2 < z < 3.5, resulting in some variation in the aver- 
age redshift across the stacked spectrum (see Figure 1). 
We ensure that the selected absorbers are not saturated 
at SDSS resolution by requiring a minimum flux of cr n /2. 
This requirement, which eliminates damped Lya systems 
and noisy lines, discards 86% of Lya absorbers in the re- 
quired optical depth range. At SDSS spectral resolution, 
solitary Lya forest lines with b— parameters ~ 30 kms -1 
and no damping wings do not reach an apparent opti- 
cal depth greater than TL yQ « 1.9, even if they are fully 
saturated. The absorption features in our sample are 
therefore blends of several Lya forest lines rather than 
individual strong lines; we address the typical column 
densities of these features in §3 below, using the stacked 
spectra themselves. 

2.2. Spectral Stacking Methodology 

For each Lya absorber in our sample, we de-redshift 
the whole spectrum to its rest-frame. Each quasar spec- 
trum is carried forward once for every Lya absorber, so 
some SDSS spectra are used more than once while oth- 
ers are not used at all. What results is a stack of spectra 
on a rest-frame wavelength grid. We employ two statis- 
tics to produce composite spectra: the median, and the 
arithmetic mean with a 3% outlier clipping. We require 
a minimum of 100 pixels for the measurement of the 
stacked spectrum at any point. When computing the 
composite spectrum redward of 1236A, we only include 
pixels if they are outside of the Lya forest of the quasar 
spectrum in question, since the pixels outside of the for- 
est have much less background absorption and associated 
random fluctuations. 

We estimate the error in the flux for the stacked spec- 
trum by bootstrapping the data using 100 realizations. 
This is mostly consistent with the estimate from the 
propagation of errors, but it is larger in the wavelength 
range between the Ly/3 and Lya lines, due to scatter in 
uncorrelated Lya forest absorption. 

Figure 1 shows the stacked spectrum we obtain using 
the median (which is our fiducial choice). The stacked 
spectrum never reaches 100% transmission (F — 1) be- 
cause there is always some uncorrelated, contaminating 
absorption. 

2.3. Composite Spectra of Lyman-a Forest Absorbers 

The full-width-half-maximum of SDSS resolution is 
136 — 167 kms -1 , which corresponds to an effective 
Doppler parameter of 6sdss = 82 — 100 kms -1 . This is 
around four times as broad as typical Lyman-series lines 
and broader still for metal lines. Hence the majority 
of contaminating "background" absorption in our spec- 
tra is not a true superposition of overlapping lines but a 
blending of distinct lines that appear superimposed due 
to SDSS resolution. 

Given that uncorrelated contaminating absorption 
varies smoothly in our stacked spectrum, we can treat 
it as a continuum by performing a standard spline fit, 
producing the 'pseudo-continuum' shown by the overlay 
in Figure 1. We have scaled away its effect by adding the 
flux decrement of the pseudo-continuum to our stacked 
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Fig. 1. — Top panel: The median stacked spectrum of 9480 
Lya absorbers. Overlaid in green is the 'pseudo-continuum' set by 
uncorrelated absorption. Lyman series lines and other lines can be 
seen along with a broad absorption trough caused by uncorrelated 
Lyman series lines. Middle panel: The mean redshift of pixels that 
make up the stacked spectrum. Bottom panel: The number of 
pixels that go into the stacked spectrum at every wavelength. 



spectrum, thus arriving at what is effectively a composite 
rest-frame spectrum of the selected Lya absorbers. 

Figure 2 shows the median spectral stack of our sample 
of Lya forest absorbers renormalized using the pseudo- 
continuum. A great many absorption lines are seen in 
this composite rest-frame spectrum. The Lyman series 
is particularly clear and is marked by dotted lines. Nu- 
merous metal lines are also seen at high confidence and 
are marked with dashed lines and labelled. 

The Lya line shows large wings, which are to be ex- 
pected as a signal of large scale structure. The scale 
of the correlation extends to ~ 3000 kms -1 , which is 
in good agreement with the findings of McDonald et al. 
(2006). There is no indication that the signal arises from 
damping wings, and higher order Lyman lines rule out 
this interpretation (sec following section). These wings 
have not been carefully fitted to separate the signal of 
clustering from the mean flux decrement in the forest, 
as the measurements in this paper are not dependent on 
this fit. 

It should be noted that the weak lines seen redward of 
some lines (e.g. O VI, Si IV, and C IV doublets) are a 
signal due to Si III interlopers. There is a small sample 
of strong Si III absorbers that has entered our sample of 
Lya lines for stacking. When these lines are misinter- 
preted as Lya, a "shadow" signal in our stacked spectra 
is seen shifted by AL yQ /Asiiii = 1216/1206. This shad- 
owing is apparent for every identified metal line in the 
arithmetic mean stack, in keeping with its status as a less 
outlier-resistant statistic. We disregard these "shadow" 
lines in our analysis below. It is surprising that no sub- 
stantial shadow is seen for Lya; we do not have a simple 
explanation for this effect, and it appears that these Si III 
absorbers are an interesting sample in their own right. 
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Fig. 2. — The rest-frame spectrum of our selected Lya absorbers. The blue vertical lines indicate the error estimate. This composite is 
produced by scaling away the 'pseudo- continuum' from the median stacked spectrum shown in Figure 1. Lyman series lines are marked 
(dotted lines). Metal lines found at high confidence are marked (dashed lines) with labels to the right. These metal lines have been fitted 
and the results of these fits are shown in Table 1. Note the difference in scales for both axes for all panels. 
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3. LINE FITTING AND INTERPRETATION 

Full interpretation of these stacked spectra will require 
the use of model spectra and/or degraded high- resolution 
spectra to account for the effects of line blending, sam- 
ple selection, and pseudo-continuum subtraction. In this 
paper, we restrict our interpretation to some general con- 
clusions that can be drawn by comparing the measured 
line strengths to simple models. 

3.1. H I Column Density 

Figure 3 shows a measurement of the H I column den- 
sity using the Lyman series. For each of seven Lyman 
series lines (n=l:Lya, n=2:Ly/3, n=3:Ly7 etc.), we mea- 
sure a rest-frame equivalent width W. We plot them 
normalized by the product /A of oscillator strength and 
line wavelength. Optically thin Lyman series lines would 
all have equal values of Wf /A. The equivalent width 
is shown for both the median and the arithmetic mean, 
with error bars drawn from bootstrapping in the compos- 
ite spectra. The noise in the arithmetic mean spectrum 
renders the pseudo-continuum fitting blueward of 935A 
unreliable, so we only show Lyman lines n=l, 2,3,4 and 
5. 

The curves in Figure 3 are created by using VPFIT 7 
to generate model lines at SDSS resolution, from which 
we measure a simulated W . The flatness of the n = 4 — 7 
points in the median stack suggests that these lines are 
optically thin, and comparison to the model curves then 
implies a total column density logiVHi/cm~ 2 w 15.3 — 
15.5 for our median absorption features. In this column 
density range, lines with velocity width b = 30kms _1 
go from significantly saturated at n — 3 to minimally 
saturated at n — 4, in agreement with the trend in the 
data points. Lines with b = 15kms , corresponding 
to the narrowest Doppler parameters seen for individual 
Lya forest absorbers, predict rising values of W/fX from 
n = 4 — 7, in clear disagreement with the observed trend. 
Lines much broader than b = 30kms _1 would predict 
too much Ly/3 and Ly7 absorption relative to the higher 
order lines. 

Of course, these absorption features are probably not 
well described by single Voigt line profiles, in part be- 
cause the composite spectrum comes from systems with 
a range of properties, but mostly because the contribut- 
ing features are themselves blends of multiple lines. The 
models in Figure 3 all under-prcdict the Lya equiva- 
lent width, plausibly because the "main" absorption fea- 
ture is near-saturated in Lya and absorption outside the 
Gaussian wings (but still within the SDSS resolution 
element) makes a larger relative contribution. While 
the b parameters implied by Figure 3 have a complex 
and non-trivial significance, the high-order Lyman series 
lines do appear to give a robust estimate of the total 
column density characteristic of our median stack, with 
logA^Hi/cm -2 fx 15.3 — 15.5. Results for the arithmetic 
mean stack suggest a column density lower by 0.1-0.2 
dex, but it is harder to draw clear conclusions because 
we cannot robustly measure Lyman-6 and Lyman-7 in 
the arithmetic mean stack. 

3.2. Metal line Column Densities 

7 http://www.ast.cam.ac.uk/ ~rfe/vpfit.html 
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Fig. 3. — Curve of growth for the first seven Lyman series lines, 
labelled 1 (Lya) to 7 (Lyrj). Points show the rest- frame equivalent 
widths normalized by /A measured from the median stack (crosses) 
and the arithmetic mean stack (squares). Curves show predic- 
tions for Gaussian lines with 6-parameters of 30 km s^ 1 (solid) and 
15kms _1 (dashed) and column densities logTVni/cm -2 = 15.1, 
15.3, 15.5 (bottom to top). 



We have used VPFIT to fit all of the metal lines 
marked in Figure 2, assuming spectral resolution R = 
2000. Wherever possible, we measure each line inde- 
pendently. Only the blended transitions C II with the 
weaker O VI line, along with O I and one of our Si II lines 
(A1304), are fitted with a joint analysis. Since contam- 
ination by Si III shadowing is a concern for the weaker 
Si IV line we discard it from the analysis. Table 1 lists 
the rest wavelengths, equivalent widths (W) and the fit- 
ted column densities (N), Doppler parameters (6), and 
reduced x values for both the median and arithmetic 
mean stacks. Many lines have fitted 6-parameters sig- 
nificantly larger than the expected width of metal lines. 
There is a large scatter in widths but low ionization lines 
are typically broader. These large velocity widths could 
be a signature of dispersion within the absorption sys- 
tems or large-scale clustering of the metal lines with the 
Lya absorbers (which is consistent with the Lya cluster- 
ing signal). These results are in line with measurements 
of C IV clustering on scales of up to 600 kms -1 (Pieri 
et al. 2006; Scannapieco et al. 2006). 

Regardless of the true line widths, these metal lines 
are likely to be optically thin, making the derived col- 
umn densities robust within the VPFIT-derived errors. 
Multiple measurements of the same species for a given 
stack are usually consistent within quoted errors, with 
the notable exceptions of Si II and C IV. In the case of 
Si II, the line at 1260A is discrepant with the other two 
measured lines. The continuum varies rapidly over a nar- 
row wavelength range in this part of the spectrum, and 
we conclude that the continuum fitting for this line is 
the dominant source of error and therefore drop the line 
from further analysis. The two lines of the C IV doublet 
imply column densities that differ by 0.1 dex, about five 
times the quoted lcr errors. Also, the fit for the strongest 
doublet member is poor. Since this is the strongest metal 
line, it may also have the strongest signal of clustering in 
the line profile, violating the single-line model assump- 
tion used in VPFIT. We conclude that the precision of 
this measurement is at the 0.1 dex level. 
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Figure 4 shows our column density measurements from 
the median and arithmetic mean stacks. Species are 
listed in order of decreasing ionization potential. The 
column densities derived from the mean stack tend to be 
slightly higher, by about 0.2 dex, but the overall agree- 
ment is good. Somewhat counter-intuitively, the simi- 
larity of line strengths in the arithmetic mean and me- 
dian stacks does not imply that this level of absorption 
is present in most selected absorbers. If the individual 
lines are not outliers with respect to the noise distribu- 
tion, then a sub-population of absorbers can shift the 
mean and median flux decrements by similar amounts. 
By looking at the full distribution of flux decrements, 
one can set limits on the fraction of systems that con- 
tribute most of the absorption. While we reserve de- 
tailed analysis to future work, our preliminary studies of 
the distribution of flux decrements at line center suggests 
that the metal lines with higher ionization potentials (for 
Si III and greater) arise in 15-30% of the selected Lya 
systems and that the lower ionization lines (for C II and 
lower) arise in 5-15% of systems. 

Despite this caveat about sample inhomogeneity, it is 
useful to compare our results with the simple predictions. 
Model curves are shown in Figure 4 using CLOUDY ver- 
sion 08.00 (Ferland et al 1998), and assuming log Am = 
15.4, a solar abundance pattern and a quasar -(-galaxy 
UV background (Haardt & Madau 2001). Four curves 
show a range of models representing typical conditions 
in the Lya forest with a metallicity of [Fe/H]= —2. 
Physical conditions p/p « 10 and T = 10 4 — 10 4 5 K 
give a reasonable match to the high-ionization lines, but 
the strengths of the lower ionization lines are under- 
predicted by two orders of magnitude or more. There- 
fore, lower ionization potential species are unlikely to be 
reproduced by models with diffuse IGM conditions. We 
can, however, get an approximate match to these lines 
with p/p ~ 1000, T w 10 4 K, and near-solar metallicities. 
Since the H I column densities implied by the Lyman- 
series lines are far too low for self-shielding (which re- 
quires Nm > 10 17 cm~ 2 ), the absorbing gas should be 
photo-ionized by the UV background, so its tempera- 
ture is unlikely to be below 10 4 K. 8 The combination 
of these high densities with Am ~ 10 154 cm~ 2 implies 
~ 10 pc pathlengths for these low ionization absorbers. 
We conclude that different absorber populations, or per- 
haps different gas phases within the same systems, are 
responsible for the high- and low- ionization lines. 

The high density, high metallicity, and short path- 
length implied by this analysis suggests that the low ion- 
ization metal lines may arise in circumgalactic regions, 
perhaps in outflows from high redshift galaxies, or per- 
haps in enriched filaments feeding these galaxies. The 
Lya absorption features are similar to those found in 
sightlines that pass within ~ 100 kpc of Lyman break 
galaxies (Steidel et al. private communication; Stcidcl 
et al. 2009). 

4. CONCLUSIONS 



8 We cannot completely rule out the possibility that the low 
ionization species arise in a small fraction of systems that are self- 
shielded (Lyman Limit systems). It is unlikely that such a model 
could be reconciled with the low mean depth of Lyman-series lines, 
but this should be investigated further. 
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log p/<p>=].0, log T = 4.5, [Fe/H] = -2 
log p/<p>=].0, log T=5.5, [Fe/H] = -2 
log p/<p> = 2.0, log T=4.5, [Fe/H] = -2 
log p/<p> = 3.0, log T=4.0, [Fe/H] = 0. 
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Fig. 4. — The column densities of metal species measured in 
order of decreasing ionization potential. Column densities from 
our median and mean composite spectrum are shown as crosses 
and squares, respectively. Models curves are shown, assuming 
logNm = 15.4, a solar abundance pattern and a quasar+galaxy 
UV background. The p/p ~ 1000 model curve is an example of 
clumpy conditions and the other curves represent densities and 
temperatures typical of the IGM as probed by the Lyo forest at 
high redshift. 



We have developed a technique for measuring the com- 
posite spectrum of Lya forest absorbers and applied it to 
13,279 quasar spectra from the SDSS. From our stacked 
spectra, we measure equivalent widths of Lyman series 
lines up to n = 7 and column densities of 19 metal lines. 
We find metal lines from species previously detected in 
the Lya forest (O VI, N V, C IV, C III, Si IV and 
Si III), from new ionization states of oxygen (O I), silicon 
(Si II), and carbon (C II), and from new elements mag- 
nesium (Mg II), aluminum (Al II) and iron (Fe II) not 
previously identified in high-redshift Lya forest studies. 
Consistent results from the median and arithmetic mean 
stacked spectra and from multiple lines of the same ionic 
species show that the measured column densities are ro- 
bust. High 6-parameters of the fitted metal lines provide 
a suggestive signal of large-scale clustering or high veloc- 
ity dispersions in the environment of the absorbers. 

Analysis of the Lyman series implies typical column 
densities Am ~ 10 15 ' 4 cm~ 2 for systems in our median 
stack. The high-ionization metal-line species can then 
be explained assuming typical physical conditions for 
the diffuse IGM and metallicity [Fe/H]« —2, but repro- 
ducing the low-ionization species requires much higher 
metallicities and higher densities, probably arising in a 
minority sub-population of the absorption systems. 

The composite spectrum technique introduced here has 
the potential to teach us a great deal about physical 
conditions and enrichment of the IGM and the spectral 
shape of the ionizing background radiation. A natural 
next step is to compare composite spectra in narrower 
bins of redshift and Lya optical depth, and to apply the 
technique to the weaker but more abundant Lya forest 
absorbers that dominate the mean opacity. In this re- 
gard, we note that the high-redshift quasar sample from 
the Baryon Oscillation Spectroscopic Survey of SDSS-III 
(Schlegel et al. 2009) will eventually exceed the size of 
the current SDSS sample by an order of magnitude. This 
technique may also be usefully applied to high-resolution 
spectra and to Hubble Space Telescope spectra of the low- 
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TABLE 1 

Metal lines measured in the composite rest-frame spectrum of Ly« absorbers 



Species 


A (A) 




Median 






Arithmetic 


mean 








logN 


o^lcms ) 


Xr (, d0t J 


W (mA) 


logN 


bfkms" 1 ) 


y'i (dof) 

/XT \ / 


W (mA) 


C II a 


1036 


13.2 ±0.4 


400 ± 300 


1.1(22) 


91 ±8 


13.3 ±0.2 


300 ± 200 


1.3 (22) 


103 ±6 


CII 


1335 


13.07 ±0.04 


260 ± 30 


1.0 (9) 


22 ±2 


13.33 ±0.02 


320 ± 20 


0.64 (9) 


38 ±2 


cm 


977 


13.13 ±0.03 


170 ± 20 


0.79 (14) 


79 ±9 


13.07 ±0.06 


250 ± 40 


1.4 (14) 


64 ±9 


CIV 


1548 


13.24 ±0.02 


160 ± 10 


5.3 (8) 


66 ± l c 


13.43 ±0.02 


180 ±10 


6.7 (8) 


100 ± 2 C 


CIV 


1551 


13.35 ±0.01 


158 ±6 


0.24 (5) 


20 ± l c 


13.53 ±0.01 


179 ±9 


1.2 (5) 


50 ± l c 


N V 


1243 


13.1 ±0.4 


150 ±20 


0.87 (17) 


9±7 


13.2 ±0.1 


210 ±80 


0.73 (17) 


14 ±8 


I b 


1302 


13.0 ±0.1 


210 ± 70 


1.0 (19) 


16 ±2 


13.49 ±0.08 


380 ± 80 


0.85 (19) 


42 ±3 


VI 


1032 


13.81 ±0.05 


210 ±30 


1.6 (7) 


76 ±5 


13.80 ±0.05 


250 ± 40 


2.9 (7) 


68 ±3 


VI a 


1038 


14.0 ±0.1 


450 ± 60 


1.1(22) 


91 ±8 


14.1 ±0.3 


550 ± 50 


1.3 (22) 


103 ±6 


Mg II 


2796 


12.3 ± 0.2 


200 ± 100 


1.7 (7) 


80 ± 20 


12.6 ± 0.1 


300 ± 100 


2.5 (7) 


150 ±20 


Mg II 


2804 


12. ±5. 


< 600 


0.98 (8) 


50 ±20 


12.7±0.3 


300 ± 100 


0.80 (8) 


70 ±20 


Al II 


1671 


11.55 ±0.04 


330 ± 40 


0.56 (11) 


14 ± 1 


11.72 ±0.03 


380 ± 40 


0.69 (11) 


20 ± 1 


Si II 


1260 


11.88 ±0.07 


200 ± 50 


0.95 (14) 


10 ±3 


12.39 ±0.06 


450 ± 90 


1.5 (14) 


37 ±3 


Si II b 


1304 


12.8 ±0.1 


230 ± 70 


1.0 (19) 


16 ±2 


13.3 ±0.1 


700 ± 200 


0.85 (19) 


42 ±3 


Si II 


1527 


12.80 ±0.05 


450 ± 60 


0.69 (16) 


16 ±2 


13.08 ±0.05 


580 ± 80 


1.4 (16) 


27 ±2 


Si IV 


1394 


12.57 ±0.01 


140 ±8 


0.52 (6) 


31 ±1 


12.75 ±0.01 


174 ±7 


0.49 (6) 


46 ± 1 


Si III 


1207 


12.66 ±0.02 


220 ± 20 


0.89 (10) 


69 ±4 


12.65 ±0.01 


253 ±8 


0.99 (10) 


66 ±3 


Fe II 


2344 


12.5 ±0.3 


500 ± 300 


1.0 (13) 


11 ± 5 


12.9 ±0.3 


< 1000 


1.3 (13) 


22 ±5 


Fe II 


2383 


12.35 ±0.09 


500 ± 100 


0.76 (19) 


37 ±7 


12.46 ±0.08 


420 ± 90 


1.3 (19) 


49 ±7 



a Lines fitted together as a blend and equivalent width measured together. 
b Lines fitted together as a blend and equivalent width measured together. 
c Partial blend so blended portion not included in equivalent width measurement. 



rcdshift Lya forest, using the power of large numbers to 
bring out features that are too weak to appear in even 
the best single-quasar spectra. 
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